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SUMMARY 

The spanwise and chordwise distribution of loading, the' lift, and 
the induced drag of a swept-back wing of slender plan form are developed 
hy means of linearized lifting-surface theory. The results are appli- 
cable for all free-stream Mach numbers. The term slender implies that 
the ratio of the reduced span (equal to the product of |l — (Mach 
number) 2 I 1 / 2 and the span) to the over-all length of the wing is small. 


INTRODUCTION 

Theoretical linearized solutions for the distribution of loading 
over lifting surfaces traveling at subsonic or supersonic speeds can be 
separated into two classes. One, in which the loading is given and the 
twisted surface required to support such a loading is found, can be 
expressed mathematically In a form necessitating the evaluation of a 
double integral involving doublets of prescribed intensity scattered 
over the plan form of the wing. The other, and usually more difficult 
problem, in which the shape of the surface is specified and the result- 
ing loading is to be determined, can be resolved into the problem of 
solving a double integral equation involving doublets of u nkn own inten- 
sity scattered over the wing plan form. It is this latter type of 
problem, sometimes referred to as a problem of the second kind, with 
which this paper is to be particularly concerned. 

In subsonic studies these double integral equations have been 
avoided by use of the lifting-line theory introduced by Prandtl. This 
simplification reduces the analysis of three-dimensional loading problems 
to the study of a single integral equation more susceptible to analysis. 
In supersonic studies, on the other hand, many important types of wing 
plan forms can be completely analyzed without further simplification. 
Outstanding examples of such plan forms are those in which the methods 
of conical flow can be used. Effectively, the presence of a conical 
flow field reduces the problem one order dimensionally 'so that again 
only a single integral equation remains to be solved. There is yet 
another type of flow pattern that will reduce the double to the single 
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integral equation arid it is proposed now to exploit this phase of the 
linearized lifting— surface problem. 

The basic linearized partial differential equation for the pej>- 
turbation velocity potential in subsonic and supersonic flow is 
known to be 


(1-M 0 2 ) $XX + $YY + ®ZZ = 0 


( 1 ) 


where Mo is the Mach number of the free stream directed parallel to 
the positive X axis. Conditions sufficient to obtain this equation from 
the more general equation based on the assumption of a perfect gas in 
nonviscous, compressible, irrotational flow are that the induced 
velocities and the nondimensional velocity gradients are small. None of 
these restrictions explicitly excludes the study of equation (1) when 
the free— stream Mach number is set equal to unity and, in fact, for wings 
with swept-back plan forms, solutions to equation (1) with Mo=l can be 
obtained that do not violate the assumptions needed to linearize the 
equation. In this case the first term of the equation vanishes. There 
are, moreover, many applications of equation (1) to aerodynamic problems 
in which the free— stream gradient of velocity $xx is itself small as 
compared to the velocity gradients in the other directions. Such is the 
usual case when the wing plan form is long and slender, that is, when 
the length of the wing is large in terms of the span. In this case, also, 
the first term of equation (1) may be neglected. 


It is possible and even preferable to combine the two conditions 
involving, in the one case, the free-stream Mach number and, in the 
other, the wing geometry into one condition that will be satisfied by 
both or either of the two. This can be done by considering the Mach 
number effect to be a stretching factor which elongates distances in the 
X direction as the Mach number ap proaches unity. In this sense a slender 

V 1 1-Mp 2 1 x span 


wing is one for which 


length 

report is valid for all such slender wings. 


is small and the theory in this 


It does not matter, then, whether M~ 1 or whether the chordwise 
gradients are comparatively small; in either case equation (1) reduces 
to 


®yy + $ZZ = 0 


( 2 ) 


vhich is simply Laplace* s equation for the perturbation velocity 
potential in a transverse plane. 


Since slender vlngs and slender ving— ' body combinations are of 
increasing practical importance, and since the assumption of slenderness 
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together with the simplicity embodied, in equation (2) permits the study 
of relatively complex aerodynamic configurations, a number of investi- 
gations based upon the resultant methods have already been published. 

In reference 1 low-aspect-ratio wings were studied and references 2, 3, 
b, 'and 5 included, respectively, a complete analysis of all the stability 
derivatives of a low^-aspect-ratio triangular wing, a lifting tri an gular 
wing with an arbitrary body of revolution, a lifting triangular cruciform 
combination on an arbitrary body of revolution, and damping— in-roll 
calculations for slender sveptr-back wings and slender wing-body combina- 
tions. In references 6 and 7 the load distribution, the lift, and the 
drag were calculated for slender swept— back wings with straight leading 
edges and with tips cut normal to the free— stream direction. The object 
of the present report is to generalize, as far as possible, the plan 
forms amenable to slender wing analysis. Thus, leading edges of rather 
arbitrary shape are included and particular attention is given to the 
effect of tip shape on the load distribution. The aerodynamic charac- 
teristics of two families of plan forms will also be included. 

A list of the important symbols used is given in the appendix. 


ANALYSIS 

The means of satisfying boundary conditions for Laplace’s equation 
are classical in nature and only a sketch of the developments to be used 
need be given here. The boundary conditions for a lifting surface 
without thickness are; 

1. Perturbation velocities vanish at infinity 

2. For all points on the Z=0 plane and not on the wing or its 
vortex sheet 

= <&z= 0 + - $Z=0— = 0 

3. For all points on the Z=0 plane 


Aft 

VdZ 


k. For all points on the Z=0 plane, within the boundaries of the 

plan form, is specified 

\ ° Z /Z=0 




= 0 


z=o- 




k 
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Applying these conditions together with the two-dimensional form of 
Green* s theorem, it is possible to express (see, for example, reference 7) 
the general solution of equation (2), in the form 1 


«(x,r,z) = — T S(X) A ?< x ,U> dY 
2n J (Y~Yi) 2 +Z' 


-S(X) 


where Y=S(X) denotes the local semispan of the wiDg. If the wing has a 
straight leading edge and if 0 is the semiapex angle, then 


S(X) = X tan 0 


In all cases the semi span of the wing will be denoted So* 

At this point in the analysis a convenient change in notation can 
be introduced. Distances in the X direction sire divided by the root 
chord C r while distances in the Y and Z directions are divided by the 
product of Cr and the tangent of the semiapex angle. Nondimens ional 
variables x,y,z are thus determined such that 


X Y 

x = — , y = , 

C r C r tan 0 


Cr tan 0 


while at the same time, by definition. 


$(X,Y,Z) = <I>(xCr, yCr tan 0, zC r tan 0) = cp(x,y, z) 


S(X) = S(xCr) = s(x) C r tan 0 


Equation (3a) is valid provided no point exists within or on the 
boundaries of the flow field such that 


lim e ^ / 0 
e — >0 ° e . 


where e is the radial distance from the point. 
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Perturbation velocities in physical space will be denoted conventionally 
by u, v , w; hence 


dx 




In the transformed variables 


- U-Cr } ^ = vCr tan 0 , ^ = wC r tan 0 

No ambiguity should arise if u, v, w are used as functional symbols 
for velocity components in either set of coordinate systems; thus 
u(X,Y,Z) transforms into u(x,y,z). 

Equation (3a) now becomes 


<p(x,y,z) 


z_ /' s(x) 

2n / ' x 
J -s(x) 


Acp( x , y ! ) dy t 

(y-yi) 2 +z 2 


(3T0) 


and this can be put in a more useful form by integrating by parts. Thus, 
since cp must be a continuous function of z everywhere except on the 
lifting surface and on the trailing vortex sheet. 


Acp(-s) = Acp(s) = 0 

and equation (3b) becomes 


Cp tan 0 p b 
2n / 

s 



( 4 ) 


where v is the perturbation velocity component in the Y direction and 
Av is the jump in the value of v in the XY plane. Since, moreover, 
s and Av are not functions of z, the derivative of both sides of 
equation (4) with respect to z leads to the result 

w ~ ~ ~ f ' 7 — — r Av dyi 
2 *J_ s (y-yi) 2 +z 2 


which as z approaches zero reduces to the final form 
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w 


z=0 



Av(x,yi)dy 1 

y-yi 


- s < y < s 


(5) 


Equation (5) is 


sometimes referred to as Prandtl's integral equation. 


Interest is now fixed on the class of problems for which the 
streamwise slope X. of the surface is given. Vertical induced velocity 
wq and X. satisfy the relation 


X(X,Y) =7T (6) 

v o 

which. simply states that the lifting surface is along a stream plane. 

The approximation has been made, of course, in equation ( 6 ) that the 
local slope of the airfoil is the same as the slope of the streamline 
passing through the Z=0 plane at the same values of. X and Y, an 
approximation usual in thin airfoil problems and consistent with the 
linearization of the partial differential equation if the slopes are 
small. Combining equation ( 6 ) with equation (5) permits the study of a 
swept-back wing of arbitrary twist and angle of attack. The case of 
greatest practical interest, however, is the one in which the lifting 
surface is a flat plate at a constant angle of attack; for this case 
w 0 is independent of X and Y. 

Solutions of the integral equation given by equation (5) are well 
known in the studies of aerodynamic problems. (See, e.g., references 8 , 
9 , and 10.) If the restriction is made that 

Av(y)dy =0 



then equation ( 5 ) can be inverted to the form 



1 

(s— y) 1 / 2 [y— (— s ) ] 1 / 2 



vo(yi)(s-yi) 1 / 2 [yi-(-s)l 1/2 

y-yi 


<tyi 


(T) 


(Equation (7) and, later, equation ( 8 ) are written in a form such that 
they can be generalized to the case where the distance to the left edge, 
— s, does not necessarily equal the distance to the right edge, s). The 
solution represented by equation ( 7 ) is extremely useful in application 
since it insures that the jump in potential at the extreme edges is 
zero. However, it is by no means the most general solution to equation 
(5). Under the restriction 
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L 


S / 


w Q (y)<iy 


__ (s-y) 1 / e [y-(-s)] 1 / 2 


= 0 


the inversion of equation (5) can also be expressed in the form 

w 0 (yi)<iyi 


Av = | {sr-y) x /*[y-{-s)] x / z J 


(y-yi)(s-y 1 ) 1 / 2 [yi-(-s)] 1 / 2 


(8a) 


to which the additional term 2 


Av = 


Aq 


(s_y)l/e[y_(_ S )]V 2 


(8b) 


can be added since this term gives zero values of Wq in the interval 
— 1 s ^ y < s for arbitrary Ao* In the solution of aerodynamic problems 
the choice of Ao is usually made by considering additional restrictions 
on the physical flow pattern such as the Kutta condition at the trailing 
edge of a lifting, surface. 

The solutions represented by equations (7) and (8) present suffi- 
ciently strong mathematical tools to provide for the study of a large 
class of swepb-back plan forms. No attempt has been made here to derive 
the most general form of the solution. It must suffice to remark that 
the inversion can be accomplished through the application of conformal 
’transformations, as was presented by Sohngen in reference 8, or through- 
a generalization of the integral operator methods commonly used in the 
inversion of Abel's integral equation. 


LOAD DISTRIBUTION FOR FLAT SWEPT-BACK WING 

The analysis derived in the preceding section is now to be applied 
to the study of swept-back wings of the type shown in figure 1. The 
family of plan forms is drawn with straight leading edges and with tips 
cut off parallel to the free-stream direction since such wings are to be 
discussed in detail. In the determination of the properties of these 


2 Terms such as 

Av = — 

(y+Ui)(s^-y 2 W 2 

(where Bi and pi are arbitrary constants (-s < pi < s)) are excluded 
by the condition given in footnote 1. 
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wings, however, it is not necessary to maintain rigid restrictions on 
the geometric parameters during the early stages of the analysis and, 
in the interest of generality, an attempt will be made to defer limita- 
tions in shape. In particular, discussion of more general types of tip 
boundaries will be included. The plan form, therefore, is assumed 
symmetrical with respect to the X axis or .free-stream direction, has a 
semiapex angle equal to 0, has a leading edge given by the equation 
Y=S(X) , has a trailing edge of unspecified shape, and is terminated 
spanwise by an arbitrary tip. 

The remainder of the treatment can best be divided into three parts 
corresponding to the three regions shown in the figure. Thus, region 1 
extends from the apex to the trailing edge of the root chord, region 2 
continues streamwise 
the remainder of the 


In region 1 the 
require w 0 to be a 
edge Acp(— s) = Acp(s) 
and the solution can 


for 

— s = y = s 


to the end of the leading edge, and region 3 includes 
wing. 

Region 1 


boundary condition! in the transformed plane 
constant for — s ^ y ^ s. Since along the leading 
= 0, equation (7) provides the complete inversion 
be written 


Av = 


2woy_ 


[ S 2(x)-y2]V 2 


(9a) 


or, if a straight leading edge is used. 


for 


Av 


gypy 

( X 2— y2)i/2 


(9b) 


= y = 


The jump in potential becomes 


Acp = Cr tan 9 


Avdyi = - 2w 0 C r tan 0 (s 2 -y2) 1 '' e 


( 10 ) 


Since the loading coefficient can be written 
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AP _ 2Au _ _2_ dA& _ 2 dAcp 
9 Vo Vo dX VoCp dx 

then 

Ap _ — Wp tan 9 s ( ds /dx ) 
q V 0 . ( S 2-y2)V 2 

or, for a straight leading edge. 


AP _ — 4vp tan 9 x 
q V 0 (xa-ya) 1 / 2 


(lit) 


Equations (11a) and. (lib) can be transformed finally to the physical 

— w Q 

plane and, setting a = — — , 


Ag _ 4aS(dS/dX) 
q (S 2 -Y 2 )V 2 


or, for a straight leading edge. 


(12a) 


AP _ l«xX tan 2 Q 
q ~ (X 2 tan 2 0-Y 2 ) 1 / 2 


(12b) 


These expressions have been derived in reference 1 for low-aspect— ratio 
wings . 


Region 2 


The boundary conditions in region 2 
the presence of the vortex sheet 
lying between the trailing edges. If, 
as in the accompanying sketch, the 
trailing edge is at +t(x), then, in 
a plane perpendicular to the x axis, 
w 0 is constant between ±t and ±s, 
and between — t and +t the value 
of A<p is consistent with the load- 
ing on upstream sections of the wing; 
that is, the strength and distribu- 
tion of vortices trailing from a 
lifting surface are proportional to ~S 
the gradient of the span loading. 


are considerably complicated by 




-t t 
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For |yj >fs|, Acp=0 . since no lifting elements occur outside the span 
of the wing. However, since a streamwise increase in Acp corresponds 
to a loading, Acp at any point behind the wing must be the same as Acp 
at the trailing edge for the same y. If the equation of the trailing 
edge is written in either of the forms 


y = t(x) 

0 

x = t*(y) 


then, for x > t* 

• Acp(x,y) •= A<p(t*.,y) 


} (13) 


(14) 


Since Acp(t*,y) is a function of the single variable y, its derivative 
with respect to y may be written in the form 


dA<p(t*,y) _ &Acp(t*,y) dt* + dAcp(t*,y) 
dy dt* dy dy 


But from the definition of 


Ap 


it follows that 


dt* 2 q J 


x=t* 


Introducing now the notation Av w for Av in the vortex wake (i.e., 
for — t < y < t) and Avp for the lateral induced velocity on the p lan 
form (i.e., for t < |y| < s), it follows that 


and 


Av p = 


1 dA<l>p(x,y) 

C r tan 9 by 


Av w 


Vq 

2 tan 9 


Ap 


q 


x=t* 


dt* 

dT + Av p 


x=t* 


(15) 


Thus the lateral induced velocity in the vortex wake is determined by 
the loading on the trailing edge and the value of Avp as the trailing 
edge is approached. 
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Using equation (5) and invoking the condition that Acp is symmet- 
rical with respect to the x axis 


_1_ P t 2yiAv v dy 1 _ _i_ p B 2y 1 Av p dy 1 
2« J q y 2 -^ 2 2n J t y 2 — y x 2 


Furthermore,, if the Kutta condition is to hold, then the loading coef- 
ficient Ap/q must he zero on the trailing edge and equation (15) 
reduces to Avy=Avpj x= t* . Thus, if t x * is written for t*(y 1; ) and 
the change in variables 


y 2 = q 

yi 2 = tjx 

2yidy x = dt)! 


Y (17) 


is made, equation (16) becomes 


_1_ T t2 Av p (t 1 *,y 1 )dT 1l 

* 2n V irm 



Avp(x,y 1 )dq I 

ni 


( 18 ) 


This is an integral equation in which w 0 is known to be a constant in 
the interval t 2 < q < s 2 . The solution to equation (18), which corre- 
sponds to an arbitrary trailing edge, will be outlined in the discussion 
of the third region. Here attention will be limited to that partic ular 
solution contained in equations (7) and (8) which will make the first 
integral of equation (18) give zero and the second integral give -2nw 0 . 
Such a solution is readily found in the. form 

Ay p * ^ (jjlq;) va f 1 ?) 


and, in subsequent developments, will prove to be of practical interest. 
Since Ap/q has been assumed zero at the trailing edge, and further, 
since Avp as given by equation (19) is zero at the trailing edge, then 
equation (15) shows that Av w must be zero across the entire vortex 
wake. Physically> such a condition implies that the span loading is 
constant when -t < y < t for when Av w is zero it follows that Aqy 
does not vary with y. This. fact will" be useful later. 
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From equation (19) the value of App can be determined as 


A<P P 

® — — 2v 0 
C r tan 6 


f 

/yi2_ t 2\ 

L 

\s2-yi 2 / 


dyi 


( 20 ) 


Introducing the notation' - 


to- - I - (l-*o 2 ) 1/2 

= dn (k Q ,u) = dh u 
s 

dy x = -sko 2 sn u cn u du 


equation (20) becomes 


Acfc, r F(*o,ko) 

—2 — = — 2w 0 sk 0 2 / cn 2 u du 

J o 


C r tan 9 


so that finally 


" Acp 


— = — 2wqS [E(\|ro>ko) — ko* 2 E(i|fo.>ko)] 


Cr tan 6 

where the elliptic integrals are defined as follows 

>0 


E(to^ko) = 


L 


(1— k 0 2 sin 2 x) l/2 dx 


E(to> k ) - f T1 2 - 2 \i/ 2 

Jo ( 1_k o 2 sin 2 x) t 


( 21 ) 


( 22 ) 


^The notation on elliptic functions is taken from Whittaker and Watson 
(reference 11) except for the use of F(i|r,k) which can be found more 
readily in tabulated form. 
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the argument \|r 0 is 


\|r 0 = arc sin 


/ s 2 — y 2 \ x / 2 

(^ S 2-t2 J 


(23) 


and the modulus ko is given by equation (21). It will be noted that, 
when t is zero, equation (22) reduces to equation (10). That is, Acp 
is continuous in passing from region 1 to region 2. This condition is, 
of course, essential. 


The loading coefficient on the wing plan form can now be found from 
equation (22) through use of the relation 


Ap _ 2 dAcp 

q VoCr dx 


In this way 


Ap 




qa tan Q 


~ **o' ~ F(\|r 0 ,k 0 ) (24) 


Again, when t is set equal to zero, the loading coefficient becomes 
that given by equation (11a). This continuity in pressure between 
regions 1 and 2 is by no means essential; in fact it will be seen later 
that in passing from region 2 to region 3 an abrupt discontinuity in 
pressure occurs. 


In deriving equation (18) the Kutta condition was assumed to hold 
at the trailing edge. This condition has yet to be applied to the 
expression for the loading coefficient in region 2 (equation (24)). So 
far, the leading and trailing edges have been expressed in arbitrary 
form. The restriction to a straight leading edge is now made and the 
trailing edge will be determined by first setting y=t and Ap/q=0 
in equation (24) and then solving for t as a function of x. This 
operation gives 


where 


dt E 0 
dx KqV 

E 0 = E^|,k 0 ^ , K 0 = F^|,k 0 ^ 


(25) 


and dt/dx is, of course, the slope of the trailing edge. The solution 
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to this equation that satisfies the condition x=l when t=0 can be 
written in the form 


t = 


Eo-ko' 2 K 0 


(26) 


That this is a solution of equation (25) can. he verified by direct 
substitution. Actually, however, the solution was obtained from 
equation (22) when the condition already mentioned that the span loading 
should be constant between the trailing edges was applied. Equation 
(26) gives t(x) explicitly as a function of the ratio t/x from which 
a graph of the trailing edge can readily be found (fig. 2). Since 
k 0 ’=t/x, equation (26) may be written in the form 


x 


1 

Eo-k 0 ,2 Ko 


(26a) 


From this form the asymptotic behavior of the trailing edge can be 
deduced, since x approaches infinity as k 0 ' approaches unity. Hence 
the trailing edge approaches asymptotically a straight line with unit 
slope, that is, becomes parallel to the leading edge. 

The asymptotic value of the chord is equal in magnitude to x-t as 
x becomes large. When equations (26) and (26a) are used, it follows 
that 


Ea-ho ,2 K 0 


which is indeterminate at k o =0. An application of L'Hospital’s rule 
gives, - however, 


(*-t) 


x— > 00 


2 

it 


This value is shown as an asymptote in figure 2. 

The solution in region 2 for the jump in potential (equation (22)) 
and the loading coefficient (equation (24)) have thus been shown to 
^apply to a wing plan form closely resembling a sweptr-back constant-chord 
‘wing except that at the root chord the trailing edge is smoothly 
filleted — a plan form of quite practical design. . (Figs. 1(a) and 1(b) 
were drawn with a trailing edge given by equation (26).) The loading 
coefficient in X,Y,Z space on a plan form with a trailing edge given by 
equation (2 6) can be written in the form 
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— = 4 tan 0 
qa 


|E(i|r 0 ,ko) 


+ 


Y 

X tan 9 


Y 2 -T 2 (X) 


I 1/2 


X 2 tan 2 0— Y 2 


_ Efi 

Ko 


F(to,ko) 



k Q = 


T 2 (X) " . 

X 2 tan 2 9 * 


to = arc sin 


” X 2 tan 2 9- Y 2 
X 2 tan 2 e-T 2 (X) 


1/2 


where T(X)=t(x)C r tan 9 is the ordinate of the trailing edge (fig. 1(a)) 
and can he determined in general from equation (26) for a given wing 
geometry and a value of X. The first pjirt of this variation for 
0 =t(x)= 1.2 is shown in figure 2. 

When the wing aspect ratio becomes very large, the loading on the 
outboard section should approach that of a two-dimensional wing, making 
an angle 9 with the X axis. Simple sweep theory (see, e.g., refer- 
ence 12) gives for the two-dimensional case the formula 


AP _ 4 tan 9 f Cn— Xn \ 

qa. (1-|3 2 tan 2 0) 1 / 2 \ X n J 


where Cn is the chord and Xn is a distance, both measured normal to 
the leading edge, and p 2 =|l-M 0 2 |. For small values of p tan 0, that 
is, under the assumptions of slender wing theory, this becomes 


®=4tane (^R Y* 
qcx \ X n J 
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It remains to compare the asymptotic behavior of equation (27) far 

outboard along the wing. To this 
end, a change is made to the coor- 
dinate system (see accompanying 
sketch) and expressed through the 



transformations 

X = Y n cos 9 + X n sin 9 + R cos 9 

Y = Y n sin 0 — X n cos 0 + R sin 0 

It is known, however, from the 
previous discussion of the asymp- 
totic value of the chord that, for 
large values of X,'T(X) approaches 
the value (X-CoJ tan 0 where 
is the value of the chord at 
infinity and k 0 approaches zero. 
Hence, as R becomes large. 


E(i)ro,0) = F(\|r o ,0) 


and in the limit 



which is in complete agreement with the result previously obtained from 
simple sweep theory. 


Region 3 

In the initial portion of the following analysis the outer boundary 
of the wing will be considered a function of X, thus 

Y = S(X) 

or 

y = s(x) 

where 

, N S(X) 

s(x) = 

C r tan 0 
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This will serve two purposes: Not only will the method he given for 
finding the solution to different tip shapes, hut also means whereby 
the trailing edge in region 2 can he modified' will he made evident. 

When the discussion is confined to region 3 > the part of the plan 
form affected hy the geometrical form of the tips (see fig. 1), it will 
he assumed that the tips are not raked in, that the trailing edge in 
region 2 is given hy equation (26), and that the leading edge in region 
2 is a straight line. Aside from the restriction just mentioned that 
raked-in tips are excluded,., the shapes of both edges in region 3 are 
for practical purposes arbitrary. An extension to the raked— in case 
will be given in the section on vortex drag at which time the variation 
of downwash in the wake is calculated. 

General solution .— The initial stages of setting up the problem are 
identical to those used in region 2. Again the Kutta. condition is 
applied at the trailing edge and equation (18) still applies. Suppose, 
now, that Avp consists of two parts such that 


Avp = Av x + Av 2 


(28) 


and set 

Avj. = [2w 0 +2g(x)] (5^) / ( 2 9) 

If Av 2 (x,y) and g(x) were both zero, this would lead to the loading 
coefficient given hy equation (24) which becomes, along the trailing 
edge y=t(x), > 


where 


and 



In this case, however, we are considering some specific plan form so that 
both t and s are fixed geometrically. Although in region 2 a differ- 
ential equation was set up for t(x) in terms of s(x) based on the 
condition that loading on the trailing edge vanishes, with the boundaries 
specified the loading in equation (30) must be considered as a residue 
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to be destroyed by the proper choice of Av 2 (x,y) and g(x). Operating 
on equation (28) to find the pressure at the trailing edge and applying 
the Kutta condition leads to the condition 


lim 
y — >t 


bf 


Y 

&V 2 (x>Yi)tyi = 2(w 0 +g) ( — Eo-ko* 


£ K °) 


ktS dt^ d* ( E ^o‘ 2 K°) 


(3D 


and Avp must, of course, satisfy equation (18) so that a further con- 
dition is supplied in the form 


_L r t2 Av 2 (t 1 *,y 1 )dT U g(x) _ _ I_ r S2 Av g (x,y^)dT ll (32) 

2n J t 0 2 71-711 2,7 J t2 


vhere t 0 , as shown in figure 1, is the distance, measured parallel to 
the y axis, from the x axis to the trailing edge in the plane which 
divides region 2 and region 3 


Equation (32) can now be inverted by the use of equation (8a), 
provided 


r 62 an r_L r t2 

J t 2 ( S 2 — T| ) - 1 / 2 ( TJ-t2 ) 1 / 2 2« J t 2 


Av 2 (t 1 *,y 1 )dT| 1 

m 


- g =0 


that is, provided 


g(x) = 




Avg(ti* > yi)dni 

( s 2 -^ 1 ) 1 / 2 ( t 2 -ri 1 ) 1 / 2 


(33) 



Applying the inversion and using this condition for g(x) 
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and 


k(ni) 


( so 2 -!!!) 1 ^ 2 


dAv2(t x *,yi) 

d^i 


(38) 



Equation (36) is a special form of Abel’s integral equation, the unique 
inversion of which is 


x Id 

h (m) = - t— 

n dr]i 



f(t2)dt2 

hx-t 2 ) 1 ^ 


as can he shown by direct substitution. 

Using the definition of h(qi) given by equation (38) and inte- 
grating with respect to t] X from t 0 2 to tj. 


Av 2 (t*,y) = - ^ Vs q 2 -t 



f (t 2 )dt 2 

( t 2 — S 2 ) J TJ-t 2 


(39) 


It is apparent that equation (33) is also a form of Abel’s in tegral 
equation and consequently can be inverted in exactly the same manner as 
equation (36). Thus, since g=0 when t=t 0 , (as Av 2 (t*,y) vanishes 
when y is zero) 


Av 2 (t*,y) 


2 



g(t 2 )dt 2 

J T) — t 2 


m 


Equaling the two values of Av 2 (t*,y) given by equations (39) and (^0) 
yields the relation 


dg_ f(t 2 ) 

dt 2 rt( t 2 — Sq 2 ) 


= 0 


which can be written, since the expression for 
equation (37) and since t 2 -so 2 = — k x 2 s 0 2 , 


2 Jig. = v o+g 
dt 2 Sq 2 — t 2 


f(t 2 ) is given by 
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This differential equation is easily solved for g ' under the condition 
mentioned above that g=0 for t=t Q so that 


g - w 0 



— tp 2 

so 2 -t 2 




Placing this in equation (4o),it follows that 


Av 2 (t*,y) 




and, if this is placed in equation (34), there results finally 


Av 2 = 2w 0 




SO 2 —!)/ 


(SQ 2 -t 0 2 )(TT-t 2 ) 

L (So 2- T|) (So 2 -^ 2 ) 


l / 2 


(41) 


The vertical induced velocity on the plan form can he now obtained 
by means of equation (28), thus • 


Avp = 2w 0 



(42) 


The velocity potential is (see the development of equation (22) from 
equation (20)) 


Aq>p 

Cr tan Q 


Avp dy x = - 2w 0 s 0 [E(iJf 2 ,k 2 ) - k 2 * 2 F(\(r 2 ,k 2 )] (43) 


so 


where 
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and, finally, since s D and t Q are independent of x, the equation for 
the loading coefficient becomes 


q 


0 


(45) 


It should be noticed that the value of Aepp is continuous with that 
given for region 2, but that the value of the loading coefficient is not. 


Discussion of Results 

The qualitative nature of the load distributions resulting from the 
application of equations (12a), ( 27 ), and (45) to a wing of triangular 
plan form and to two wings of swept— back plan forms is shown in figure 3. 
The loading in the case of the triangular wing is, of course, the same as 
that gi'ven in reference 1 and, in fact, differs from the exact solution of 
the linearized equation for arbitrary Mach number only by t he factor l/ E* 
where E’ is the complete elliptic integral with modulus J 1— (3 2 tan 2 9 . 

The deviation in this case can therefore be assessed accurately. The 
results in figures 3( b ) and 3(c) can be compared in a more general manner 
with the load distributions for a constant-chord wing obtained in reference 
13 for supersonic Mach numbers. These results, shown in figure 4, are 
qualitatively quite similar to those obtained by means of slenden-wing 
theory. Across the Mach cone from the root— chord trailing edge the 
loading in both cases is continuous, and falls rapidly in the after portion 
until it reaches zero at the trailing edge of the wing. Across the Mach 
cone from the leading edge of the tip the loading is discontinuous and 
behind it the magnitude of the loading is close to zero. The principal 
difference between the results of the present paper and those of refers 
ence ,13 appears in the position of the discontinuities, which, in slender^ 
wing theory, occur at the Mach lines corresponding to a free— stream Mach 
number of unity, rather than at Mach lines determined from exact free— 
stream conditions. 


AERODYNAMIC CHARACTERISTICS OF FLAT SWEPT-BACK WINGS 

Span Loading 

The span loading for the wings shown in figure 3 can best be studied 
in terms of the circulation function T(y). In magnitude the circulation 
equals the jump in the velocity potential at the trailing edge of the wing 

r(Y) = (46) 


and total lift L is expressed in the form 
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rtSo pT.E. riS o 

L = / dY / ApdX = pV 0 / T(Y)dY (47) 

S 0 "'L.E. "'-So 

where L.E. and T.E. stand for leading and trailing edges. Figure 5 
shows the spanwise variation of P for the three types of plan forms 
shown in figure 3> each having the same wing root chord but with types 
(b) and (c) having larger spans. From equation (10) circulation for the 
type (a) wing is 

T(Y) = 2V 0 a ( S 0 2 — Y 2 ) 1 /2 (48) 


while for the two other cases equation (22) gives, for 0 ^ |Y| < 

To. 


and, for T 0 g |Y | g S Q , 


r(Y) = 2V 0 aS 0 


r(Y) = 2VoaSo[E(^ 2 ,k 2 ) - k« 2 2 F(fe,k 2 ) ] 


(^9a) 


(4 9 b) 


where 




arc sin 


( 


Sq2~Y2 N 

Sq 2 — Tq 2 / 


1/2 

) 


k* 2 


= (l-k 2 2 ) 1/2 


Tfi 

So 


The shapes of the curves in figure 5 reveal the differences in the 
basic characteristics of the three wings. 'First, the triangular wing has 
an elliptic span loading and therefore will have the least induced drag 
for a given span of the three wings. Second, the type (b) wing has a 
span loading close to the elliptic and should have characteristics quite 
similar to the triangular wing. Third, wing (c) has the same span loading 
as wing, (b), but at the same time has considerably more wing area. The 
constant span loading given in equation (49a) corresponds, of course, to 
the physical conditions implicit in the functional form of Avp as 
introduced in equation (19)* 


Lift 

The lift on the three types of wings can be obtained from direct 
integration. Substituting from equations (48) and (49) into equation 
(47)’ it follows that 



NACA TH 1992 


25 


L = - 4voPVo(T 0 Cr tan 0 + Ij.) 


(50) 


where I x is the double integral given by 



/ 


pEE 

of So 2 -' 


- Y 2 
■To 2 


k P 2 



d| 


(51) 


But I x can be evaluated by making the substitution Y=|i S Q and 
changing the order of integration. Thus 


Ii = S Q 2k 2 2 


/> (ghri 


,(i-k 2 2 i 2 ) i/2 


<^1 


Integrating with respect to and using the Jacobian transformation 

|=sn(u,k 2 ) 

Ii = So 2 k 2 2 ^ — So 2 k* 2 (E 2 -^£ , 2 2 K 2 ) 


Finally, since the function r given by equations (48) and (4-9) is 
continuous, then at Y=T 0 , 


So(®* 2 — k* 2 2 K 2 ) — Cr tan 0 


and 


Ii — So 2 k 2 2 ^ — ToCp tan 0 (52) 

This gives for lift (equation (50)) 


| = 2rtoS 0 2 



(53) 


If the. results for lift are to be given in coefficient form, the 
area of the wing plan forms must be calculated. Figure 6, obtained after 
a numerical integration, can be used to obtain wing aspect ratio as a 
function of semispan for given root chord and semiapex angle. The aspect 
ratio A is, of course, given by 
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A = ( 2S °) 2 
S 


where S is wing area. From equation (53 ) > lift coefficient is 
expressible in the form 


CL _ « A 
a tan 0 2 tan 0 



(54a) 


Equation -(54a) is valid for wings of types (b) and (c) and includes 
type (a) as a special case. With the information given by figures 2 
and 6, the value' of Cl/<x tan 0 can be plotted against A/tan 0 for 
these wings as in figure 7. It should be stressed that equation (54a) 
holds also when T 0 is zero and reduces in that case to 


(Cl)t o =0 = | Aa (5 1 *) 

9 

When To equals zero for the type (b) wing, the plan form is triangular 
and equation (54b) agrees with the result for that case in reference 1. 
When T o =0 for the type (c) wing, the trailing edge of the root chord 
is behind the leading edge of the tip and, as shown in figure 7, lift 
coefficient is then a linear function of aspect ratio. This transition 
in the lift coefficient occurs at A=2.48 tan 0. ' 

When the aspect ratio of either the type (b) or (c) wing becomes 
large, it approaches the value 4S 0 2 tan 0/(So 2 — To 2 ) and the lift coeffi- 
cient becomes , 


(Cl) A ->oo = 2na tan e 


(54c) 


Equation (54c) corresponds to the exact expressions for lift coefficient 
of a swept— back inf inite— aspect— ratio wing at a free— stream Mach number 
equal to 1 and also agrees to the first order with the low— speed lift 
coefficient of a highly swept, infinite— aspect— ratio wing. 

Vortex Drag 

The formula for the vortex drag is well known from the study of 
incompressible-flow theory where it is referred to as the induced drag. 

Its derivation for both subsonic and supersonic theory depends on the 
calculation of the momentum transport through a plane perpendicular to the 
X axis and infinitely far behind the wing. 



NACA TN 1992 


27 


The vortex drag of the wings under consideration can be calculated 
from the expression 


D V = - 


pCr tan 9 


/° 

J -So 


(v 0 ) TdY 


(55) 


where (w 0 ) is v 0 at X=» and where, from equation (5), 


-So 


(56) 


From equation (49a), for 0 ^ |Yi| ^ To 


df 

dYn 


= 0 


(57a) 


and from equation (49b), for T 0 = |Yi| i Sq 


JLT 

dY 


r /Yi^o 2 V /2 Yi 


(57b) 


Placing these expressions in equation (56) , rearranging, and using the 
transformation Yi = So*f\Ly there results 


( \ V ° a 


(To/So) 172 


dTji 

T-Tll 


Hi— ( T o/Sp) 

1— T|1 


2 V2 


(58) 


and, after direct integration. 


— V 0 a 


2 ,,2 v 1/2 


\So 2 -Y 2 J 


, 0 < Jy| < To 


< v oL = < 


/ (59) 


-JJ 0 a , T 0 < |Y | g S 0 


V 
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This result is shown in the accompanying sketch for types (h) and (c) 

wings and for triangular wings, the 

latter being included as a special 

case in the analysis when To=0. 

Since the induced velocities in an 

X = constant plane are a function 

only of the vortices streaming through 

it (equation (5)), the value of w 

behind the wing is the same in every 

plane. The quantity (w 0 ) , there— 

00 

fore, represents the' downwash in 
regions immediately behind the wing 
as well as at points infinitely 
distant. With this concept and the 
results presented in the sketch, we 
are immediately able to verify tne 
loading results for type (c) wings, 
the loading over type (b) wings being 
given, since the sketch shows that 
.the downwash behind the portion of 
the trailing edge, which is cut normal 
to the free stream, is the same as the 
downwash. on the wing itself. Hence, 

the boundary condition that w=w 0 i n the plane of the wing is automatically 
satisfied for any variation in plan form directly behind this edge, and no 
loading is required to force the streamlines into a pattern which they 
already follow. 



The effect of the instability and resultant rolling up of the vortex 
sheet is such that the magnitude and distribution. of downwash. can be used 
in computing induced effects off the wing only a short distance back of 
the trailing edge. However, since the total kinetic energy is not 
changed by this rolling up, the drag of the wing can be determined from 
the velocities computed on the basis of the undistorted vortex sheet in 
the XY plane. Hence, vortex drag can now be determined from equation 
(55) with T as given in equations (49a) and (49b) and with (wo)^ 
obtained from equation (59)* The expression for vortex drag then becomes 


D v = 2pV 0 2 a 2 C r 



dY + 


2pV 0 2 a2i 1 


where I x has already been determined (equation (52)). The elliptic 
integral can be evaluated by the transformation Y=So dn (u,k 2 ) where, 
as in equation (49b), 


k» 2 


To 

Sq 
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The final expression for drag is 

~ = S 0 2 k 2 2 « - 4s 0 C r tan 0 (E«2-k 2 2 K» 2 ) 


(60) 


where K * 2 and E * 2 are complete elliptic integrals of the first and 
second kind with modulus k’ 2 . In coefficient form, equation (60) becomes 


CPy _ A ^ k 2 2 it _ E , 2 -k 2 2 K t 2 \ 
a 2 tan 0 tan 0 4 . So/Cp tan 0 y v 

A plot of equation ( 6 l) is shown in figure 8 for all three types of 
wings. For large aspect ratios the vortex drag decreases as aspect ratio 
increases and, in fact, it can be shown that Cpy approaches zero as 
aspect ratio becomes infinitely large. However, it is also apparent from 
the figure that the factor Cd v /(Cl 2 /«A) increases slightly with aspect 
ratio in the range shown. It is the latter factor that is minimized by 
the elliptic span loading of the triangular wing. 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Moffett Field, Calif., Oct. 18 , 1949. 


APPENDIX 

TABLE OF IMPORTANT SYMBOLS 


A 

Cr 

Cl 

Cd v 

D v 

EqjEi 


aspect ratio of wing 

root chord of wing 

lift coefficient of wing (L/qS) 

vortex drag coefficient of wing (D v /qS) 

vortex drag of wing 

complete elliptic integrals of the second kind with moduli k 0 ,k 1 , 
respectively 

incomplete elliptic integral 'of the second kind 

. „ 0.1/2 

(1— k sin x) dx 



E(*,k) 
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incomplete elliptic integral of the first kind 
ilr 

‘ T dx 

J (1— k 2 sin 2 x) 1 / 2 _ 


function introduced in equation (29) 


(-8) 


_ If. \ 1/2 

S2 ) 


_ tf_v 

so 2 y 




M) 


1 - ^ 


T rt \l/2 


(l_ k 2)V2 

complete elliptic integrals of the first kind with moduli kpjkj, 
respectively 

lift of wing 

free— stream Mach number ; 


static pressure 

1 

free— stream dynamic pressure 


(i° v ° 2 ) 


y coordinate of wing tip or wing leading edge (See fig. 1(h).) 
[s(x) C r tan 9 = S(X) ] 

Y coordinate of wing tip or wing leading edge (See. fig. 1(a).) 

particular value of s(x) determined by wing— tip location 
(See fig. 1(b).) ( s 0 C r tan 9 = S 0 ) 

particular value of S(X) determined by wing-tip location 
(See fig. 1(a).) 

wing area 

y coordinate of wing trailing edge (See fig. 1(b).) 

[t(x) Cr tan 9 = T(X) ] 


Y coordinate of wing trailing edge (See fig. 1(a).) 
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to 

To 

u(x,y,z) 

v(x,y,z) 

v(x,y,z) 

u(X,Y,Z) 

v(X,Y,Z) 

v(X,Y,Z) 

V 0 

( w o) a> 

x,y,z 

X,Y,Z 

a 

3 

r 

q 

0 

p 


particular value of t(x) determined by wing-tip location 
(See fig. 1(b).) (t c C r tan 0 = T 0 ) 

particular value of T(X) determined by wing— tip location 
(See fig. 1(a).) 


transformed values of perturbation velocities in X, Y,Z direc- 
tions 


perturbation velocities in X,Y,Z directions 

free— stream velocity 

vertical induced velocity at X=oo, Z=0 

transformed Cartesian coordinates 

(xC r = X, yC r tan 0 = Y, zC r tan 0 = Z) 

Cartesian coordinates in physical space 

angle of attack of wing 

*/ 1 1-Mo 2 | 

circulation 

y 2 

semiapex angle of wing 
free-stream density 


cp(x,y,z) transformed perturbation velocity potential 
[$>(xC r , yC r tan Q, zC r tan 0)] 

$(X,Y,Z) perturbation velocity potential 
Ap 

loading coefficient 

pressure on lower surface — pressure on upper surface 

q 

Av lateral velocity discontinuities at Z=0 

1 

Av x , A v 2 components of Avp 

Acp discontinuity in perturbation potential at z=0 


; ) 
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Subscripts 

p pertaining to wing plan form 

v pertaining to vortex wake 
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Figure 2 Plot showing trailing edge as given by equation (26). 





Figure 5- Span wise variation of section circulation for the wing plan 
forms of figure 3. 






Figure 6.- Variation of aspect ratio with wing span for swept-back wings having trailing edges 
given by equation (26). 



a tan 0 





tan 0 

Figure 7- Variation of lift-curve slope with aspect ratio. 









